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ABSTRACT 

This  paper  is  focused  on  the  study  of  effect  of  heat  and 
mass  transfer  on  chemically  reacting  boundary  layer  flow 
of  a  Casson  fluid  over  a  porous  stretching  sheet  in  the 
differential  equations  are  transformed  by  introducing 
similarity  variables  and  solved  numerically  by  using 
Range  Kutta  method  along  with  shooting  method.  The 
velocity,  temperature  and  concentration  distributions  are 
discussed  number  tables  for  different  parameters  entering 
into  the  problem.  It  is  observed  that  the  Casson 
parameter  decreases  the  velocity  field  while  the 
temperature  is  enhanced  with  increasing  Casson 
parameter. 

Keywords:  Mass  Transfer,  MHD,  Casson  Parameter, 
suction,  heat  generation,  chemical  reaction,  radiation 
effect  and  Soret  effect 

Introduction: 

The  combined  effects  of  heat  and  mass  transfer  with 
chemical  reaction  have  attracted  many  researchers  due  to 
its  wide  range  of  applications  in  Engineering  and 
Science.  Heat  and  mass  transfer  occur  simultaneously  in 
the  processes  of  drying  and  evaporation  at  the  surface  of 
a  water  body,  energy  transfer  in  a  wet  cooling  tower  and 
the  flow  in  a  desert  cooler.  The  study  of 
magnetohydrodynamic  (MHD)  flow  of  an  electrically 
conducting  fluid  is  of  considerable  interest  due  to  its 
application  in  many  engineering  problems  such  as  MHD 
generators,  plasma  studies,  nuclear  reactors,  and 
geothermal  energy  extractions.  Das  et  al.  (1994)  studied 
the  first  order  chemical  reaction  effect  on  the  flow  past 
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an  impulsively  started  infinite  vertical  plate  with 
constant  heat  flux  and  mass  transfer.  Anjalidevi  and 
Kandaswamy  (1999)  considered  the  heat  and  mass 
transfer  on  steady  laminar  flow  along  a  semi-infinite 
horizontal  plate  in  the  presence  of  chemical  reaction. 
The  Soret  and  Dufour  effects  on  heat  and  mass 
transfer  about  vertical  surfaces  in  porous  media  have 
been  studied  by  Postelnicu  (2004).  radiation  on 
boundary  layer  flow  over  a  moving  vertical  porous 
plate  was  analyzed  by  Makinde  (2005). 
Muthucumaraswamy  et  al.  (2006)  have  studied  the 
effect  of  chemical  reaction  isothermal  vertical  plate 
with  radiation.  Mandal  (1985)investigated  mass 
transfer  effects  on  unsteady  hydromagnetic  free 
convection  flow  past  an  accelerated  vertical  porous 
plate.  Jha  (1991)  studied  the  effect  of  magnetic  field 
on  mass  transfer  flow  past  a  uniformly  accelerated 
vertical  plate  through  a  porous  medium.  Elbashbeshy 
(1997)  analyzed  the  heat  and  mass  transfer  along  a 
vertical  plate  in  the  presence  of  magnetic  field.  The 
combined  heat  and  mass  transfer  convection  flow 
from  a  vertical  surface  with  Ohmic  heating  and 
viscous  dissipation  was  analyzed  by  Chen  (2004). 

The  study  of  boundary  layer  flow  over  a  stretching 
sheet  has  finds  applications  in  chemical  engineering, 
particularly  in  manufacture  artificial  film,  artificial 
fibers,  polymer  extrusion,  drawing  of  plastic  films  and 
wires,  glass  fiber  and  paper  production.  Crane  (1970) 
was  the  first  researcher  who  investigated  the  boundary 
layer  flow  over  a  stretching  surface.  After  the  study  of 
fluid  flow  over  a  stretching  sheet  has  received  wide 
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attention  among  researchers.  Rajagopal  et  al.  (1984) 
discussed  the  flow  of  second  order  fluid  over  a  stretched 
sheet.  Anderson  et  al.  (1992)  consider  magnetic  field  on 
the  flow  of  a  viscoelastic  fluid  past  a  stretching  sheet. 
Abel  et  al.  (2005)  analyzed  MHD  boundary  layer  flow 
over  continuously  moving  stretching  surface  embedded 
in  a  porous  medium  by  considering  the  Buoyancy  force 
and  effects.  Mukhopadyaya  et  al.  (2008)  discussed  the 
free  convective  boundary  layer  flow  with  variable 
viscosity  over  a  stretching  surface  with  thermal  radiation. 
Pal  (2009)  investigated  the  mixed  convection  flow  of  an 
incompressible  fluid  over  a  stretching  sheet  in  the 
presence  of  radiation.  Ahmed  (2009)  analyzed  the  free 
convective  heat  and  mass  transfer  of  a  viscous 
incompressible  fluid  over  a  stretching  sheet  in  the 
presence  of  suction  with  Soret  and  Dufour  effects.  Heat 
and  mass  transfer  in  non-Newtonian  fluids  have 
applications  in  engineering  such  as  catalytic  reactors,  the 
filtration  and  blood  plasmapheresis  devices.  Casson  fluid 
is  the  most  popular  non-Newtonian  fluid  used  to  model 
blood.  It  is  defined  as  a  shear  thinning  liquid  which  has  a 
infinite  viscosity  at  zero  rate  of  shear,  possessing  a  yield 
stress  below  which  no  flow  occurs  and  zero  viscosity  at 
infinite  rate  of  shear.  It  is  reduced  to  Newtonian  fluid  at 
very  high  wall  shear  stress  i.e.  when  the  wall  stress  is 
much  greater  than  yield  stress.  Merill  et  al.  (1965)  and 
McDonald  (1974)  conducted  experiments  on  the 
behavior  of  blood  as  a  Casson  fluid.  Eldabe  (1995) 
considered  the  heat  transfer  of  Casson  fluid  flow 
between  two  rotating  cylinders.  The  flow  of  Casson  fluid 
in  a  tube  was  studied  by  Dash  et  al.  (2000)  and  Nagarani 
et  al.  (2004).  Mass  transfer  in  a  Casson  flowing  through 
an  annular  geometry  was  examined  by  Nagarani  et  al. 
(2006).  Attia  (2010)  analyzed  the  transient  Couette  flow 
of  a  Casson  fluid  between  parallel  plates  with  magnetic 
field  and  heat  transfer.  The  unsteady  boundary  layer  flow 
of  a  Casson  fluid  over  a  moving  flat  plate  was  studied  by 
Mustafa  et  al.  (2011).  Hayat  et  al.  (2012)  studied  the 
mixed  convection  stagnation  point  flow  of  a  Casson 
fluid.  Shehzad  (2013)  discussed  the  effects  of  mass 
transfer  on  the  MHD  boundary  layer  flow  of  a  Casson 
fluid  with  chemical  reaction. 

In  view  of  the  above  studies,  we  consider  the  MHD 
boundary  layer  flow  of  a  Casson  fluid  over  a  stretching 
sheet  with  heat  and  mass  transfer  and  chemical  reaction. 

.Nomenclature: 

u,v  :  Velocity  components  along  the  x  and  y 

axes 

T  :  Fluid  temperature  inside  the  boundary 


layer 

T 

1  co  • 

the  plate 

Ambient  temperature  for  away  from 

T 

W 

wall 

Uniform  constant  temperature  at  the 

C  :  Species  concentration  inside  the 

boundary  layer 

fluid 

Species  concentration  of  the  ambient 

g  : 

Acceleration  due  to  gravity 

k  ; 

Permeability  of  the  porous  medium 

Bo  : 

Uniform  magnetic  field 

Pr  : 

Prandtl  number 

M  : 

Magnetic  parameter 

Gr  : 

Grashof  number 

Gc  : 

Modified  Grashof  number 

Sc  : 

Schmidtnumber 

D  : 

concentration 

Molecular  diffusivity  of  the  species 

K  : 

Permeability  parameter 

Greek  Symbols: 

^ :  Dimensionless  temperature 

P  :  Fluid  density 

^ :  Kinematics  viscosity 

:  Electrical  conductivity 

n  0  • 

? 

coefficient 

Thermal  and  concentration  expansion 

a  . 

Thermal  diffusivity 

F 

w 

Suction  parameter 
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Mathematical  Formulation: 


Concentration  Equation: 


Consider  the  steady,  incompressible  flow  of  a  Casson 
fluid  over  a  porous  stretching  surface  at  y  =  0.  Choose 
the  coordinate  system  such  that  x-axis  is  parallel  to  the 
surface  and  y-axis  normal  to  the  surface.  The  fluid 
occupies  half  space  y  >  0.  A  uniform  magnetic  field  BO 
is  applied  in  the  y  direction.  The  transverse  applied 
magnetic  field  and  magnetic  Reynolds  number  are 
assumed  to  be  very  small,  so  that  the  induced  magnetic 
field  and  Hall  effects  are  negligible.  We  also  considered 
the  heat  and  mass  transfer  processes  in  the  presence  of 
chemical  reaction  and  suction.  The  rheological  equation 
of  state  for  an  isotropic  and  incompressible  flow  of  a 
Casson  fluid  can  be  written  as  (Nakamura  and  Sawada 
(1988),  Mustafa  et  al.  (2012)). 

_ 2(Pb  +  py  1 V2 x)ey  ,7r  >  nc 
lJ  2 (juB  +  Py  /  4ln)eij ,  n  <  nc 


dc  dc  _  _  s2c 

- b  V —  D  - 

U  dx  dy  8y2 


d  T 

-MC-cj  +  a  — 


<5/ 


(4) 


^-  =  4az(T-TJ 

dy 


Where  u  and  v  are  the  velocity  components  in  x  and  y 
directions,  y  is  the  kinematic  viscosity, 


P  Mb  V  71  c  y  js  qie  Casson  fluid  parameter, a  is 
the  electric  Conductivity  of  the  fluid, p  is  the  density 
of  the  fluid,  T  is  the  temperature  of  the  fluid,  C  is  the 
Concentration  field,  a  is  the  thermal  diffusivity,  D  is 


the  mass  diffusivity,  k  i  is  the  reaction  rate. 


The  Boundary  Conditions  for  the  Velocity, 
temperature  and  Concentration  fields  are 


T=T  C=C 

u=u  w,v=- V 0  5  w  when  y=0 


Where  pB  is  the  plastic  dynamic  viscosity  of  the  non- 
Newtonian  fluid,  Py  is  the  yield  stress  of  the  fluid,  e  ij 
denotes  the  (i,  j)  -  th  component  of  the  deformation  rate, 
71  =  e  ij  e  ij  is  the  product  of  the  component  of 
deformation  rate  with  itself,  7ic  is  the  critical  value  of  ji 
based  on  the  non-Newtonian  model. 


Under  above  conditions  the  continuity,  momentum  and 
energy  equations  are 

Continuity  Equation: 


du  dv 
— + —  =  0 
dx  dy 


(1) 


u-mw,  T->T  »,C->C“  as  y->  oo  (5) 

We  define  the  following  similarity  variables  as 


\-yJin ) 

y  x 


¥ 


T,  -  T, 


Ain) 


Cw-Ca 


B  = 


CpCf 


(6) 


Where  ¥  is  the  stream  function  with  u= 

dy/  _  dy/ 

dy  ’  dx 


Momentum  Equation: 


du  du  1  d2u 

u - bv —  =  r(lH — ) — ^ —  <t 

dx  dy  P  dy2 

Bl  y 

— u-—ru  +  gfi(P~Poo)  +  g  rfrr  r  \ 
P  K  p  {Ls-Lsn) 


Energy  Equation: 


dT  dT  d  T 

u - 1- v — -a — — 

dx  dy  dy2 


1  dq^ 

pC  P  dy 


+  Q^(T-T) 
PC,  ^ 


(2) 


(3) 


^is  the  similarity  variable  substituting  the  non- 
dimensional  variables 

In  view  of  (2),  (3)  and  (4)  the  Equations  (2)  to  (4)  take 
the  form 


+r -if')2  - (M + K)f  =  -{GtP + Gap) 

(7) 

—0"  +  f0,-(Ra-B)0  =  0 

Pr  V  '  (8) 
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f  -  Scf  9'  -  ScK  ,<j)  +  9"SrSc  =  0  ^ 

where  the  primes  denote  the  differentiation  with  respect 
tor),  ^is  the  magnetic  parameter,  K  is  the  permeability 
parameter,  ^'is  the  temperature  Grashof  number,  is 

the  concentration  Grashof  number,  Ptis  the  Prandtl 
number,  Ra  is  the  Radiation  parameter,  B  is  the  Heat 

Generation,  Sr  is  the  Soret  effect  and  &  is  the  Schmidt 
number. 

The  Corresponding  non  -dimensional  boundary 
conditions  are 

/  =  U  /  =  S,0  =  \,<t>  =  1,  when^=0 

/'(°o)->  O,0(oo)->  0  ^(oc)  _»  o  asy->oc  (10) 

Solution  of  the  Problem 

The  governing  boundary  layer  equations  (7)  to  (9) 
subject  to  boundary  conditions  (10)  are  solved 
numerically  by  using  shooting  method.  First  of  all, 
higher  order  non-linear  differential  equations  (7)  to  (9) 
are  converted  into  simultaneous  linear  differential 
equations  of  first  order  and  they  are  further  transformed 
into  initial  value  problem  by  applying  the  Range  Kutta 
method  along  with  shooting  technique.  From  the  process 
of  numerical  computation,  the  skin-friction  coefficient, 
the  Nusselt  number  and  Sherwood  number  which  are 

respectively  proportional  ft)  ^),  -  ^(O)  ancj  -^(0) 
are  also  sorted  out  and  their  numerical  values  are 
presented  in  a  tabular  form. 

Results  and  Discussion 

In  order  to  get  a  physical  insight  into  the  problem,  a 
representative  set  of  numerical  results  is  shown 
graphically  in  Figs.  1-27,  to  illustrate  the  influence  of 
physical  parameters  viz.,  the  effect  of  the  Casson 
parameter  (P),  Magnetic  field  (M),  Porous  parameter  (K), 
Grashof  Number  (Gr),  Modified  Grashof  number  (Gc), 
Suction  parameter  (S),  Prandtl  Number  (Pr),  Radiation 
parameter  (R),  Heat  Generation  (B),  Schmidt’s  Number 
(Sc),  Chemical  Reaction  (Kl)  and  Soret  effect  (Sr)  on 
the  flow  variables  have  been  graphically  presented. 

Fig  1  presents  the  effect  of  Casson  parameter  (P)  on 
velocity.  It  is  observed  that  the  velocity  increases 
asymptotically  from  its  highest  value  on  the  surfaces  to 
zero  as  r|— »go.  The  presence  of  yield  stress  reduces  the 
velocity.  Increasing  values  of  p  increases  the  velocity 
further  and  thus  there  is  a  decrease  in  the  thickness  of  the 


boundary  layer. 

Fig.2  shows  the  variation  of  chemical  reaction  (Ch) 
on  velocity.  The  effect  of  chemical  reaction  increases 
with  the  increases  on  velocity.  Fig  3  illustrates  the 
effect  of  Grashof  number  (Gr)  on  velocity.  We 
observe  that  the  velocity  decreases  with  the  increase 
of  Grashof  number  (Gr)  with  take  the  reverse  action. 
Fig  4  illustrates  the  effect  of  modified  Grashof 
number  (Gc)  on  velocity.  We  observe  that  the  velocity 
decreases  with  the  increase  of  modified  Grashof 
number  (Gc) 

with  take  the  reverse  action.  Fig  5  illustrates  the  effect 
of  porous  parameter  (K)  on  velocity.  We  observe  that 
the  velocity  decreases  with  the  increase  of  porous 
parameter  (K)  with  take  the  reverse  action.  Fig  6 
illustrates  the  effect  of  magnetic  field  (M)  on  velocity. 
We  observe  that  the  velocity  decreases  with  the 
increase  of  magnetic  field  (M)  with  take  the  reverse 
action.  Fig  7  illustrates  the  effect  of  radiation 
parameter  (R)  on  velocity.  We  observe  that  the 
velocity  decreases  with  the  increase  of  radiation 
parameter  (R).  Fig  8  shows  the  temperature  for 

different  values  of  heat  generation  (B).  It  is  observe 
that  the  temperature  increases  with  the  increase  of 
heat  generation  (B).  From  fig  9  shows  the  temperature 
for  different  values  of  Casson  parameter(P).  We 
observe  that  the  temperature  increases  with  the 
increase  of  Casson  parameter  (P).  Fig  10  shows  the 
temperature  for  different  values  of  chemical  reaction 
parameter  (Ch).  It  is  seen  that  the  temperature 

decreases  with  the  increase  of  chemical  reaction 
parameter  (Ch).Fig  11  shows  the  temperature  for 
different  values  of  modified  Grashof  number  (Gc).  It 
is  seen  that  the  temperature  decreases  with  the 
increase  of  modified  Grashof  number  (Gc).  Fig  12 
shows  the  temperature  for  different  values  of  Grashof 
number  (Gr).  It  is  seen  that  the  temperature  increases 
with  the  increase  of  Grashof  number  (Gc).  Fig  13 
shows  the  temperature  for  different  values  of  porous 
parameter  (K).  It  is  seen  that  the  temperature 

decreases  with  the  increase  of  porous  parameter  (K). 
Fig  15  shows  the  temperature  for  different  values  of 
Prandtl  number  (Pr).  It  is  seen  that  the  temperature 
decreases  with  the  increase  of  modified  Prandtl 
number  (Pr).Fig  15  shows  the  temperature  for 

different  values  of  radiation  parameter  (Ra).  It  is  seen 
that  the  temperature  decreases  with  the  increase  of 
radiation  parameter  (Ra).  Fig  16  shows  the 

concentration  for  different  values  of  heat  generation 
(B).  It  is  seen  that  the  temperature  decreases  with  the 
increase  of  heat  generation  (B).Fig  17  shows  the 
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concentration  for  different  values  of  Casson  parameter 
(P).  It  is  seen  that  the  concentration  increases  with  the 
increase  of  Casson  parameter  (P).Fig  18  shows  the 
concentration  for  different  values  of  chemical  reaction 
(Ch).  It  is  seen  that  the  concentration  decreases  with  the 
increase  of  chemical  reaction  (Gc).Fig  19  shows  the 
concentration  for  different  values  of  modified  Grashof 
number  (Gc).  It  is  seen  that  the  temperature  decreases 
with  the  increase  of  modified  Grashof  number  (Gc).Fig 

20  shows  the  concentration  for  different  values  of 
Grashof  number  (Gr).  It  is  seen  that  the  concentration 
increases  with  the  increase  of  Grashof  number  (Gr).Fig 

2 1  shows  the  concentration  for  different  values  of  porous 


parameter  (K).  It  is  seen  that  the 
concentrationincreases  with  the  increase  of  porous 
parameter  (K).Fig  22  shows  the  concentration  for 
different  values  of  magnetic  parameter  (M).  It  is  seen 
that  the  concentration  increases  with  the  increase  of 
magnetic  parameter  (M).Fig  23  shows  the 
concentration  for  different  values  of  Schmitt  number 
(Sc).  It  is  seen  that  the  concentrationdecreases  with 
the  increase  of  Schmitt  number  (Sc).Fig  24  shows  the 
concentration  for  different  values  of  Soret  number 
(Sr).  It  is  seen  that  the  concentration  decreases  with 
the  increase  of  Soret  number  (Sr). 


Fig  1:  The  Velocity  Profile  for  different  values  for  Fig  2:  The  velocity  profile  for  the  different  values 
Casson  Parameter  (P)  for  chemical  reaction  (Ch) 


Fig  3:  The  velocity  profile  for  the  different  values  for 
modified  Grashof  number  (Gc) 


Fig  4:  The  velocity  profile  for  the  different  values 
for  Grashof  number  (Gr) 
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Fig  5:  The  velocity  profile  for  the  different  values  Fig  6:  Xhe  velocity  profile  for  the  different  values 
for  porous  parameter  (K)  for  magnetic  field  (M) 


Fig  7:  The  velocity  profile  for  the  different  values  Fig  8:  The  temperature  profile  for  the  different 

for  radiation  parameter  (R)  values  for  heat  generation  (B) 


Fig  9:  The  temperature  profile  for  the  different 
values  for  Casson  parameter  (P) 


T1 


Fig  10:  The  temperature  profile  for  the  different 
values  for  Chemical  Reaction  (Ch) 
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Fig  11:  The  temperature  profile  for  the  different  Fig  12:  The  temperature  profile  for  the  different 
values  for  modified  Grashof  number  (Gc)  values  for  Grashof  number  (Gr) 


Fig  13:  The  temperature  profile  for  the  different  Fig  14:  The  temperature  profile  for  the  different 
values  for  porous  parameter  (K)  values  for  Prandtl  number  (Pr) 


Fig  15:  The  temperature  profile  for  the  different 
values  for  radiation  parameter  (Ra) 


Fig  16:  The  Concentration  profile  for  the  different 
values  for  heat  generation  (B) 
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Fig  17:  The  Concentration  profile  for  the  different 
values  for  Casson  parameter  (P) 


1.5 

h 


3 


Fig  18:  The  Concentration  profile  for  the  different 
values  for  Chemical  reaction  (Ch) 


Fig  19:  The  Concentration  profile  for  the  different 
values  for  modified  Grashof  number  (Gc) 


Fig  20:  The  Concentration  profile  for  the  different 
values  for  Grashof  number  (Gr) 


Fig  21:  The  Concentration  profile  for  the  different 
values  for  porous  parameter  (K) 


Fig  22:  The  Concentration  profile  for  the  different 
values  for  magnetic  parameter  (M) 
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Fig  23:  The  Concentration  profile  for  the  different 
values  for  Schmitt  number  (Sc) 

CONCLUSIONS: 

This  paper  gives  the  effect  of  Heat  and  mass  transfer  on 
MHD  boundary  layer  flow  for  a  Casson  fluid  over  a 
stretching  sheet  in  presence  of  heat  generation  and 
chemical  reaction.  The  governing  equations  are  solved 
by  shooting  technique.  It  is  observed  that  velocity  is 
decreasing  function  of  Casson  parameter  p,  Magnetic 
parameter  M,  Suction  parameter  S.  The  thermal 
boundary  layer  thickness  decreases  with  an  increase  in 
Prandtl  number  Pr.  The  effect  of  P  and  M  on  Skin 
friction  coefficient  is  opposite. 

>  We  observe  that  the  velocity  decreases  with  the 
increase  of  Grashof  number  (Gr)  with  take  the 
reverse  action. 

>  We  observe  that  the  velocity  decreases  with  the 
increase  of  modified  Grashof  number  (Gc)  with  take 
the  reverse  action. 

>  It  is  observe  that  the  temperature  increases  with  the 
increase  of  heat  generation  (B). 

>  We  observe  that  the  temperature  increases  with  the 
increase  of  Casson  parameter  (P). 

>  It  is  seen  that  the  temperature  decreases  with  the 
increase  of  heat  generation  (B). 

>  It  is  seen  that  the  temperature  decreases  with  the 
increase  of  modified  Grashof  number  (Gc). 

>  It  is  seen  that  the  concentration  increases  with  the 
increase  of  Grashof  number  (Gr). 

>  It  is  seen  that  the  concentration  decreases  with  the 
increase  of  Soret  number  (Sr). 


Fig  24:  The  Concentration  profile  for  the  different 
values  for  Soret  effect  (Sr) 
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